Universal critical exponent in class D superconductors 
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We study a physical system consisting of non-interacting quasiparticles in disordered supercon- 
ductors that have neither time-reversal nor spin-rotation invariance. This system belongs to class D 
within the recent classification scheme of random matrix ensembles (RME) and its phase diagram 
contains three different phases: metallic and two distinct localized phases with different quan- 
tized thermal Hall conductances. We find that critical exponents describing different transitions 
(insulator-to-insulator and insulator-to-metal) are identical within the error of numerical calcula- 
tions and also find that critical disorder of the insulator-to- metal transition is energy independent. 

PACS numbers: 73.20.Fz, 72.15.Rn 



The properties of quasiparticles in disordered super- 
conductors have been intensively investigated. Those 
systems are representatives of new symmetry classes dif- 
ferent from the familiar three classes in normal disordered 
conductors and in the Wigner-Dyson random matrix en- 
sembles. A list of additional random matrix ensembles, 
determined by these new symmetry classes has been es- 
tablished about 10 years agoi. These additional ran- 
dom matrix ensembles describe zero-dimensional prob- 
lems, and are appropriate to model a small grain of a 
superconductor in the ergodic highly conducting limit. 
In the corresponding higher-dimensional systems of the 
same symmetry classes there can be transitions be- 
tween metallic, localized, or quantized Hall phases for 
quasiparticles^^. 

The symmetry class D (which we address in this paper) 
may be realized in superconductors with broken time- 
reversal invariance, and either broken spin-rotation in- 
variance (as in d-wave superconductors with spin-orbit 
scattering) or spinless or spin-polarized fermions (as in 
certain p-wave states). The associated changes in quasi- 
particle dynamics must be probed by energy transport, 
since neither charge density nor spin are conserved. In 
this paper we present extensive numerical results on a 
model first introduced by Cho and Fisher (CF)S (see a 
detailed description below) , which has a particularly rich 
phase diagram in two dimensions^. We use advanced 
numerical calculations, proposed in Ref. [5| and described 
in detail in&£ to overcome round-off errors in calcula- 
tions of renormalized localization lengths. We also apply, 
for the first time, an optimization algorithm to deter- 
mine both critical exponent and critical disorder for the 
insulator-to-metal transition and critical exponent for the 
insulator-to-insulator transition. Our results suggest the 
existence of the universal critical exponent v = 1.4 ± 0.2 
for the thermodynamical localization length £ for both 
types of transitions mentioned above, an energy inde- 
pendent value of the critical disorder W cr « 0.19 ± 0.02 
and show collapse of all data on one curve with thermo- 
dynamic localization length £ ~ [e(W cr — W / )]~ l/ . 

The original network modet^ was proposed to describe 
transitions between plateaux in the quantum Hall effect 
(QHE). QHE is realized in a two-dimensional electron 



gas subjected to a strong perpendicular magnetic field 
and a random potential. When random potential varies 
smoothly (its correlation length is much larger than the 
magnetic length) a semiclasscial description becomes rel- 
evant: electrons move along the lines of constant poten- 
tial. When two equipotential lines come close to each 
other (near a saddle point), tunneling is feasible. In the 
network model, electrons move along unidirectional links 
forming closed loops in analogy with semiclassical mo- 
tion on contours of constant potential. Scattering be- 
tween links is allowed at nodes in order to map tunnel- 
ing through saddle point potentials. Propagation along 
links yields a random phase <fi, thus links are presented 
by diagonal matrices with elements in the form exp(i</>). 
Transfer matrix for one node relates a pair of incoming 
and outgoing amplitudes on the left to a corresponding 
pair on the right; it has the form 



cosh sinh 9 
sinh cosh 9 



(1) 



The node parameter < 
in the following way 



is related to the electron energy 



■ ln(sinh 9), 



(2) 



where e is a relative distance between the electron en- 
ergy and the barrier height. It is easy to see that the 
most "quantum" case (equal probabilities to scatter to 
the left and to the right) is at e = (9 — 0.8814), in fact, 
numerical calculations show that there is an extended 
state at that energy. 

Numerical simulations on the network model are per- 
formed on a system with fixed width M and periodic 
boundary conditions in the transverse direction. Multi- 
plying transfer matrices for N slices and then diagonal- 
izing the resulting total transfer matrix, it is possible to 
extract the smallest Lyapunov exponent A (the eigenval- 
ues of the transfer matrix are exp(AiV)). The localization 
length £m is proportional to 1/A. Repeating calculations 
for different system widths and different energies it is 
possible to show that the localization length satisfies 
a scaling relation 



6/ = f (M_\ 



(3) 
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When renormalized localization length £m /M becomes 
M-independent it strongly suggests that the system is 
in the critical state. Indeed, it means that when you 
double the system size M, the localization length £m 
also doubles, or, from the mathematical standpoint, the 
only way to make Eq. (3) M-independent and consistent 
is to put the thermodynamic localization legth £ equal to 
infinity (extended state). 

For a D class symmetry a Bogoliubov - de Gennes 
Hamiltonian may be written in terms of a Hermitian 
matrb*i. The corresponding time evolution operator is 
real, restricting the generalized phase factors to be O(N) 
matrices for a model in which TV-component fermions 
propagate on links, and to the values ±1 for N = 1, the 
case that was studied. There are three models: random 
bond Ising mode l 11 ' 12 , supporting two different localized 
phases, uncorrelated 0(1) model^, where phases on the 
links are independent random variables and all states are 
extended 5 , and the model first introduced by Cho and 
Fisher (CF)£ where scattering phases with the value tt 
appear in correlated pairs. Each model has two param- 
eters: the first one is a disorder concentration W, such 
that there is a probability W (1 — W) to have a phase 
(tt) on a given link. The second parameter is an energy 
e describing scattering at the nodes. For the CF model, 
the phase diagram (updated version of which is presented 
in Fig. 3) in the e-W plane has three distinctive phases: 
metallic, and two insulating phases characterized by dif- 
ferent Hall conductances. The sensitivity to the disorder 
is a distinctive feature of class D. 

The existence of region of extended states means that 
the smallest Lyapunov exponent at each particular en- 
ergy is zero or extremely small. This point was discussed 
in detail in Ref. [j| The advanced algorithm using the 
structure of the transfer matix imposed by current con- 
servation and the symmetry of class D suggested there 
allowed to increase the accuracy of the data significantly. 
However, in the vicinity of the critical points the noise 
increases and in order to identify the values of critical 
disorder W cr at different energies one needs to focus only 
on large system sizes M , therefore errors in the values 
of localization lengths £m can still be significant. In this 
paper we have decided to address this problem from the 
other direction. We use the values of we can trust 
(not very close to the critical point) and demand that 
all renormalized localization lengths satisfy scaling rela- 
tions (allowing the program to find the optimal values of 
critical exponent and critical disorder). 

In the CF model the disorder is introduced only at 
the nodes, allowing for the offdiagonal elements of Eq. 
(1) to be multiplied by ±1 (disorder probablity W is the 
probability of that factor to be —1). We have started 
numerical calculation for fixed values of W changing the 
system widths M and the energies e. In those cases an 
obvious critical energy is e — (the most " quantum" case 
explained above) and we have been looking for the criti- 
cal exponent v which fits all the data onto one curve. We 
use a special optimization program which checks different 
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FIG. 1: (a) Renormalized localization lengths for various 
system widths M and energies e for a fixed disorder value 
W = 0.1. (b) Scaling of the same data with critical exponent 
v = 1.4 



critical exponents and chooses the optimal one. The typ- 
ical data for the renormalized localization lengths for dif- 
ferent system widths M and energies e for a fixed disorder 
FK = 0.1 is shown in Fig. la. For large e the states are lo- 
calizaed almost ideally (renormalized localization length 
is proportional to 1/M). As one approaches the critical 
energy e = 0, the values of all renormalized lengths in- 
crease significantly (tending to infinity). In Fig. lb it is 
shown that all the data fit onto one curve when we scale 
the horizontal axis as eM 1 /" with the critical exponent of 
the thermodynamic localization length v = 1.4. Indeed, 
substituting into the right-hand side of Eq. (3) £ ~ <l~ v 
one can get the variable of the horizontal axis in Fig. lb. 
We have performed such calculations for different fixed 
values of the disorder W in the range [0.05; 0.15]. The 
optimization program have found very close values of v 
for all the data. 

Next, we address the most interesting problem - defin- 
tion of the critical disorder and corresponding critical 
exponent. We have fixed the value of energy e and vary 
M and W . The renormalized localization lengths are 
shown in Fig. 2a for a fixed e = 0.1. We then use the 
optimization program to find optimal values for both the 
critical disorder W cr and critical exponent. The results of 
the scaling are shown in Fig. 2b. Surprisingly, the value 
of the critical exponent found in this way was identical 
(within the error of numerical calculations) as for a fixed 
disorder. The critical disorder found for a wide range 
of different fixed values of energy e (e = 0.1, 0.2, ...1.) 
has the same value W cr independent on e. This way of 
defintion of W cr by the otimization procedure looks more 
promising, especially taking into account that when one 
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FIG. 2: (a) Renormalized localization lengths for various 
system widths M and disorder W for a fixed energy value 
e = 0.1. (b) Scaling of the same data with critical exponent 
v = 1.4 and critical disorder W cr = 0.19. 

looks for £m/M to become M-independent in the stan- 
dard procedure, the errors increase significantly and one 
needs to take M as big as possible and to trust only to 
the largest values of M. We, therefore, suggest that our 
method leads to the phase digram shown in Fig. 3 with 
the vertical critical line separating insulator and metal 
(for energies e not very close to zero). We, however, can- 
not determine the exact position of the multicritical point 
at e = (studied previouslj^^) , which is underlined by 
the dots in Fig. 3. 




FIG. 3: Updated phase diagram for a CF model with metallic, 
insulating and quantized Hall phases. 

The very fact that we have obtained almost identi- 
cal values of the critical exponent v = 1.4 ± 0.2 for 
both types of transition, insulator-to-insulator transi- 
tion with £ ~ (T v and insulator-to-metal transition with 
£ ~ (yV cr ~W)~ v (which is rather different from anything 
obtained before, (e.g. two different critical exponents for 
a class CSii ) has led us to the last fit. We have suggested 



FIG. 4: Collapse of all of the data on a single curve with uni- 
versal critical exponent v = 1.4 and universal critical disoredr 
W cr = 0.19. 



that the thermodynamic localization length should de- 
pend on the product of both distances of energy and dis- 
order from the critical ones with the same critical expo- 
nent in the following way £ ~ [e(W C r — VF)] - ". Indeed, 
all the data for the renormalized localization lengths ob- 
tained for various e, W and M collapse on the same curve 
after suggested fit as shown in Fig. 4. 

Detailed analytical explanations for the divergence of 
the localization length in various types of class D sys- 
tems for a one-dimensional model (M = 2) were pre- 
sented earlier— Here we wish to add a qualitative ar- 
gument for a two-dimensional system. First, consider a 
system without disorder (W = 0). Then, there are no 
phases on the links, tunneling alone defines a quaiparti- 
cle behaviour, and for any non-zero energy makes system 
an insulator. For any e ^ at any node quasiparti- 
cle prefers to turn into the same direction circling closed 
trajectories. Futhermore, because the probability of a 
tunneling event is proportional to exp(— 7re), it is trivial 
to derive that after N tunneling events the probability 
becomes exp(— vreiV), and, considering a number of tun- 
neling events TV as a distance, one immediately obtains 
the localization length as £ ~ e . We have found that 
those statements, mentioned above, are correct numeri- 
cally as well. For W = and e ^ all states are ideally 
localized (£m is M-independent) and fit onto one curve 
depending on £/M with £ ~ e . Now let us introduce 
disorder. It causes the appearence of phase 7t on some 
links. Those phases can cause a constructive interfer- 
ence for different trajectories leading from one point of 
the network model to another, making therefore an ex- 
tended state. For a 0(1) model phases on the different 
links are uncorrelated, and therefore even a small disor- 
der W introduces that constructive interference immedi- 
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ately. In the Cho-Fisher model that we consider in this 
paper, phases on both sides of the node are correlated 
(the same), therefore small disorder does not mean an 
immediate dclocalization, but it is clear that the stronger 
the disorder is the larger is the opprotunity for the state 
to be extended. 

To summarize, we have studied critical exponents de- 
scribing the divergence of thermodynamic localization 
lengths for two both types of transition (insulator-to- 
insulator and insulator-to-metal) of the CF model and 
found that within the error of numerical calculations both 
exponents are equal to v = 1.4 ± 0.2. We have used op- 
timization procedure in order to determine not only the 
best fit of the data leading to the critical exponent, but 
also to look for a critical disorder if the energy is fixed. 



We have found that the value of the critical disorder is 
energy independent (for energies e not very close to zero). 
Obviously, we cannot rule out completely some possible 
weak dependence of the critical disorder on the energy, 
especially when e — > oo, but for a wide range of energies 
studied, the critical disorder is W cr — 0.19 ± 0.02. 
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